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(N 

^ ■ rations of quantum gravity on the propagation of "flavoured" (Klein-Gordon and Dirac ) neutral 



In the first part of this work we discuss possible effects of stochastic space-time foam configu- 



OO . particles, such as neutral mesons and neutrinos. The formalism is not the usually assumed Lind- 



which the propagation of the matter particles is considered. We arrive at expressions for the re- 
spective oscillation probabilities between flavours which are quite distinct from the ones pertaining 



\ blad one, but it is based on random averages of quantum fluctuations of space time metrics over 

00 

o 

\^ \ to Lindblad-type decoherence, including in addition to the (expected) Gaussian decay with time, 

o ■ 

a modification to oscillation behaviour, as well as a power-law cutoff of the time-profile of the 
respective probability. In the second part we consider space-time foam configurations of quantum- 
D , fluctuating charged black holes as a way of generating (parts of) neutrino mass differences, mim- 

icking appropriately the celebrated MSW effects of neutrinos in stochastically fluctuating random 



^ ■ media. We pay particular attention to disentangling genuine quantum-gravity effects from ordinary 

effects due to the propagation of a neutrino through ordinary matter. Our results are of interest 
to precision tests of quantum gravity models using neutrinos as probes. 

PACS numbers: 03.65.Yz, 14.60.Pq, 14.60.St, 04.70.Dy 



1 



I. INTRODUCTION AND MOTIVATION 



The important feature of classical General Relativity, is the fact that space-time is not 
simply a frame of coordinates on which events take place, but is itself a dynamical entity. For 
conventional quantisation this poses a problem, since the space-time coordinates themselves 
appear "fuzzy". The "fuzzyness" of space-time is associated with microscopic quantum 
fluctuations of the metric fleld, which may be singular. For instance, one may have Planck 
size (10"'^^ m) black holes, emerging from the quantum gravity (QG) "vacuum", which may 
give space-time a "foamy" , topologically non-trivial structure. 

An important issue arises which concerns the existence of a well-defined scattering matrix 
in the presence of black holes, especially such microscopic ones (i.e for strong gravity); the 
information encoded in matter fields may not be delivered intact to asymptotic observers. 
In this context we refer the reader to a recent claim by S. Hawking according to which 
information is not lost in the black hole case, but is entangled in a holographic way with the 
portion of space-time outside the horizon. It is claimed that this can be understood formally 
within a Euclidean space-time path integral formulation of QG. In this formulation the 
path-integral over the topologically trivial metrics is unitary, but the path integral over the 
topologically non-trivial black hole metrics, leads to correlation functions that decay to zero 
for asymptotically long times. Consequently only the contributions over trivial topologies 
are important asymptotically , and so information is preserved. In simple terms, according 
to Hawking himself, the information is not lost but may be so mangled that it cannot be 
easily extracted by an asymptotic observer. He drew the analogy to information encrypted 
in "a burnt out encyclopedia", where the information is radiated away in the environment, 
but there is no paradox, despite the fact that it is impossibly difficult to recover. 

However, there are fundamental issues we consider as unanswered by the above interesting 
arguments. This makes the situation associated with the issue of unitarity of effective 
matter theories in foamy space-times unresolved. On the technical side, one issue that 
causes concern is the Euclidean formulation of QG. According to Hawking this is the only 
sensible way to perform the path integral over geometries. However, given the uncertainties 
in analytic continuation, it may be problematic. Additionally, it has been argued l| that 
the dynamics of formation and evaporation of (microscopic) black holes is unitary using 
Maldacena's holographic conjecture of AdS/CFT correspondence j2| for the case of anti-de- 
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Sitter (supersymmetric) space-times. This framework describes the process in a very specific 
category of foam, and may not be vahd generally for theories of QG. However even in this 
context the role of the different topological configurations is actually important, a point 
recently emphasised by Einhorn 0]. In Maldacena's treatment of black holes P], the non- 
vanishing of the contributions to the correlation functions due to the topologically non-trivial 
configurations is required by unitarity. Although such contributions vanish in semiclassical 
approximations, the situation may be different in the full quantum theory, where the role 
of stretched and fuzzy (fluctuating) horizons may be important, as pointed out by Barbon 
and Rabinovici M. 

n 

The information paradox is acutest in the case of gravitational collapse to a black hole 
from a pure quantum mechanical state, without a horizon; the subsequent evaporation due 
to the celebrated Hawking-radiation process, leaves an apparently "thermal" state. It is in 
this sense that the analogy is made with the encoding of information in the radiation 
of a burning encyclopedia. However the mangled form of information in the burnt out 
encyclopedia, is precisely the result of an interaction of the encyclopedia with a heat bath 
that burned its pages, thereby leading to an irreversible process. The information cannot 
be retrieved due to entropy production in the process. 

In our view, if microscopic black holes, or other defects forming space-time foam, exist in 
the vacuum state of quantum gravity (QG), this state will constitute an "environment" which 
will be characterised by some entanglement entropy, due to its interaction with low-energy 



matter. This approach has been followed by the authors jo, in many phenomenological 

n 



tests or microscopic models of space-time foam gj], within the framework of non-critical 
striug theory; the latter, m our opm.ou. is a viable (,ron-equihbriu,n) theory of spaee-fme 

foam '9], based on an identiflcation of time with the Liouville mode. The latter is viewed as 
a dynamical local renormalization-group scale on the world-sheet of a non-conformal string. 
The non-conformality of the string is the result of its interaction with backgrounds which 
are out of equilibrium, such as those provided by twinkling microscopic black holes in the 
foam. The entropy in this case can be identifled with the world-sheet conformal anomaly of 
a o"-model describing the propagation of a matter string in this fluctuating background Qj. 
Although within critical string theory, arguments have been given that entanglement entropy 
can characterise the number of microstates of Anti-de-Sitter black holes Q], we do not flnd 
these to be entirely convincing. 
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In view of the above issues, it is evident that the debate concerning space-time foam 
remains open. The thermal aspects of an evaporating black hole are suggestive that the 
environment due to quantum-gravity is a sort of "thermal" heat bath. This has been pursued 
by some authors, notably in ref. Another proposal, the D-particle foam model j^, 

considers the gravitational fluctuations that could yield a foamy structure of space-time 
to be D-particles (point-like stringy defects) interacting with closed strings. There are no 
thermal aspects but there is still the formation of horizons and entanglement entropy within 
a fluctuating metric framework. 

In general, for phenomenological purposes, the important feature of such situations is the 
fact that gravitational environments, arising from space-time foam or some other, possibly 
semi-classical feature of QG, can still be described by non-unitary evolutions of density 
matrices. Such equations have the form 

9tp = Aip + A2P (1.1) 

where 

AiP = ^ [p, H] 

and H is the hamiltonian with a stochastic element in a classical metric. Such effects may 

n 

arise from back-reaction of matter within a quantum theory of gravity which decoheres 
the gravitational state to give a stochastic ensemble description. Furthermore within mod- 
els of D-particle foam arguments in favour of a stochastic metric have been given . The 
Liouvillian term A2P gives rise to a non-unitary evolution. A common approach to A2P, 
not based on microscopic physics, is to parametrise the Liouvillian in a so called Lindblad 
form We note at this point that any non-linear evolutions that may characterise a 

full theory of QG (see e.g. a manifestation in Liouville strings ^^), can be ignored to a first 
approximation appropriate for the accuracy of contemporary experimental probes of QG. 
Generically space-time foam and the back-reaction of matter on the gravitational metric 
may be modelled as a randomly fluctuating environment; formalisms for open quantum me- 
chanical systems propagating in such random media can thus be applied and lead to concrete 
experimental predictions. The approach to these questions have to be phenomenological to 
some degree since QG is not sufficiently developed at a non-perturbative level. 

One of the most sensiti ve p robes of such stochastic quantum-gravity phenomena are 



neutrinos 
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I9I |20| , in particular high-energy ones 



2l\ . It is the point of this 



article to present various approaches to gravitationally-induced decoherence of matter and 
to classify some characteristic experimental predictions that could be falsified in current or 
near future neutrino experiments. 

The neutrino, being almost massless, and weakly interacting, can travel long distances 
in the Universe essentially undisturbed. Thus the detection of high energy neutrinos, which 
are produced at early stages of our Universe, say in Gamma-Ray-Bursters or other violent 
phenomena, can carry important information on the Universe's past which would not have 
reached us otherwise. If space-time has therefore a stochastic foamy structure, the longer the 
neutrino travels the greater the cumulative quantum-gravity effects become. For instance, 
due to their known mass differences, the neutrinos exhibit oscillations between their various 
flavours, and such oscillations appear to attenuate with time in stochastic environments. 
Although such an attenuation may be too small to be detected in laboratory experiments, 
it may nevertheless be appreciable in the case of ultra-high-energy neutrinos, which have 
travelled cosmological distances before reaching the observation point on Earth [2l,|2l|. From 
such (non) observations of damping effects, one may place important bounds on quantum- 
gravity effects, information that may prove quite useful in our theoretical quest of under- 
standing space-time. 

Moreover, there is another interesting possibility regarding neutrinos. As pointed out 
recently in llG*, the tiny mass differences between neutrino flavours may themselves (in 
part) be the result of a CPT violating quantum-gravity background. The phenomenon, if 
3e the generalisation of the celebrated Mikheyev-Smirnov-Wolfenstein (MSW) 
The latter arises from effective mass differences between the various neutrino 
flavours, as a result of different type of interactions of the various flavours with matter within 
the context of the Standard Model. The phenomenon has been generalised to randomly 



true, wou 



effect m 



fluctuating media 



24| , which are of relevance to solar and nuclear reactor /^-decays neutrinos. 



This stochastic MSW effect will be more relevant for us, since we consider space-time foam, 
as a random medium which induces flavour-sensitive mass differences. 

The structure of the article will be the following: we commence our analysis by consid- 
ering in sec. II flavour oscillations between two generations of neutrinos, whose dynamics 
are governed by Klein-Gordon or Dirac Lagrangians in the presence of weakly fluctuating 
background random gravitational fields. The Klein-Gordon case is an idealisation when the 
effects of neutrino spin are ignored. Moreover it can be of interest in its own right when 
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flavour oscillations of neutral mesons are considered. The case of Dirac particles with two 
flavours is considered in section III. An effective description in terms of two-level systems is 
derived and analysed. We then proceed in sec. IV to discuss gravitational MSW effects in 
oscillation phenomena (also for two flavours) for the case when the particles are highly rel- 
ativistic (a situation applicable to neutrinos). We pay particular attention to disentangling 
potential genuine quantum-gravity-induced decoherence effects from conventional effects due 
to the passage of the neutrino probe through ordinary stochastic fluctuating matter. As we 
shall discuss, the disentanglement is achieved via the energy E and oscillation length L 
dependence of the relevant probability. In particular, conventional effects attenuate to zero 
as the parameter L/E ^ 
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26l |. in contrast to the genuine quantum-gravity decoherence 



effects which, at least in some models of space-time foam decoherence, exhibit a L ■ E de- 
pendence. Conclusions and outlook are presented in section V, followed by three appendices 
that contain some technical details of our formalism. 



II. GRAVITATIONAL DECOHERENCE CALCULATIONS FOR SCALAR PAR- 
TICLES 

Since the effects of stochastic space-time foam can appear through both Aip and A2P 
in (jU we shall for clarity isolate their individual signatures. The most satisfactory way of 
dealing with the effects of such a background is by coupling covariantly the gravitational fleld 
to a Klein-Gordon or Dirac lagrangian.This avoids intuitive arguments which are sometimes 
presented j27| and correctly incorporates covariance unlike these other approaches. 

For the case of scalar particles of mass m, such as neutral mesons (or in the toy case where 
the spin of a neutrino of mass m is ignored), we can describe the motion of the particle in a 
curved background by means of a Klein-Gordon equation for a fleld $. The Klein-Gordon 
equation in a gravitational fleld reads: 

g'^'^D^Dp^ - = 0. (2.1) 

where the metric tensor and is a covariant derivative. We will consider the neutrino 
to be moving in the x-direction. For simplicity [27I1 we will examine the situation where 
the relevant part of the contravariant metric can be regarded as being in 1 + 1 dimension. 
Moreover if metric fluctuations are caused by D-particle foam ^ there are further arguments 



in favour of such a truncated theory. A small stochastic perturbation of the flat metric can 
be written as 

g = 07]0^ (2.2) 

with 

O^f-^' - V r,J-''] (2.3) 

and where the static coefficients Oj's are gaussian random variables satisfying (oj) = and 
{aittj) = 5ij(Ji. This is a simplified model and could be made more complicated, for example, 
by having a general symmetric covariance matrix for the a^'s. Such complications will not 
affect our qualitative results and magnitudes of estimates. From ()2.2|) : 

^^u^l -{ai + lf + al -a3(ai + l) + 02(04 + 1) I 4) 
I -03(01 + 1) + 02(04 + 1) -o| + (04 + 1)^ / 

Since the Christoffel symbols FJ^^ = and i? = for static Oj's the Klein-Gordon equation 
is 

{g'^^dl + 2g''d,d, + g''dl)<P - = 0. (2.5) 
For positive energy plane wave solutions 

0(x,t) ~ ^(A;,«;)e*(-"*+^") 

we have the dispersion relation 



,01 

UJ 



+ Zl^y^id'^n^ - g^'ig^'k^ + M2). (2.6) 
For an initial a flavour state with momentum k, the density matrix p at time t is 



Pit) = E Uo.jU;^U:iU,,e^(-'-^^V^){f,\. (2.7) 

where /3 is a flavour index and j, l{= 1, 2) denote indices for mass eigenstates with eigenvalue 
M = mi and M = m2.The bras and kets in 12. 71 are flavour eigenstates (corresponding to the 
flavours denoted by the subscripts) and U is the mixing matrix which can be parametrised 
by an angle 6: 

cos 9 sin6' , 

U=\ I (2.8) 

— sin 6 cos 6 



Now since the tu's are functions of classical random variables (which thus have a positive 
probability distribution), the averaging of p{t) over these random variables is a positively 
weighted (generalised) sum over density matrices. Hence the averaged density matrix is also 
positive and represents a mixed state. The probability of transition from an initial state of 
flavour 1 to 2 is 

Prob(l 2) = ^ t/yf/*j.f/*;f/2ie*('^'-^^^* (2.9) 



where the time dependent part is 



Since the {at} are assumed to be independent Gaussian variables, our covariance matrix S 
has the diagonal form 

/ ^ \ 
0^00 
0^0 

V 

with ffj > 0. The calculation of transition probabilities requires the evaluation 



(2.10) 



)= / (i^aexp(— a ■ S ■ a)e 



i(uji—u>2)t 



detS 



From ()2.fij) we obtain 
1 



UJl — 002 
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00 



(^01)2p _ ^00(^11^2 + ^2) _ ./(^01)2p _ ^00(^lip + ^2) (2.12) 



(2.11) 



Now, since fluctuations are small, we can make the expansion 

_1_(^(^01)2^2_ ^00(^11^2 + ^2)) ^ ^ (^^) + ^ (^^) ^. + ^ J.^. (^^) + 



(2.13) 
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where the non-zero expansion coefficients are 

c{mi) 



= -^A;2 + mf 



di{mi) = -^/WTrr^, d^imi) 



fc2 



fu{mi) = yp+T^, /mK) = -^^=^ 



(2.14) 



f22{mi) 



mf+2k'^ 



, f2i{rni) 



JUymi) — 2 (ifc2+„i2)3/2 



and /j, is symmetric. In this approximation we find that 



where 



i{(jJl—UJ2)t\ 



- exp f — I exp(i6t). 

VX2/ 



B 



i6i 








0-2 





-i(<i-*' 

-ifc2St 

2d 







2d 
J_ 








^ - \ik^~ct j 

(T4 2 / 



(2.15) 



Pi = + 2id6 (^/e^ - d) + b^k^aicr^t^, 
P2 = 4(i2 - 2i£ {li^ZaA + 26ai) t + O (a^) 



with 



6 = ypT^- x/F" 



m 



21 



c = ml{k'^ + ml) — m^A;^ + ml) 



d= ^k'^ + ml^W^ 



rrio. 



It is particularly illuminating to consider the limit k » mi, 777.2 for which d — k^, b — 
where (Am)^ — ml — ml, and c = ^^^^ . We then have 



(2.16) 



(Am)2 
2k ' 
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PlP2 = 



) 




) 



Hence we see that for highly energetic scalar particles the stochastic model of space-time 
foam leads to a modification of oscillation behavior quite distinct from that of the Lindlbad 
formulation. In particular for the transition probability there is a gaussian decay with 
time, a modification of the oscillation period as well an additional powerlaw fall-off both 
decays are invariant under t ^ —t which is of course related to their origin from Ai. From 
this characteristic time dependence bounds can be obtained for the fluctuation strength of 
space-time foam. They are compatible with previous estimates and will be discussed later. 

III. DECOHERENCE OF DIRAC PARTICLES 

Although scalar flavour oscillation is the relevant case for neutral mesons, for the im- 
portant case of neutrino oscillations and space-time foam it can only be a rudimentary 
approximation. The spinorial structure should be incorporated into the description. The 
usual discrete level descriptions of oscillation phenomena cannot suggest the natural way 
to incorporate the background and this leads to consideration of the Dirac equation in the 
presence of a stochastic gravitational background. For definiteness we will take neutrinos 
to be described by two flavours and by massive Dirac spinors \E'; also a term is introduced 
which incorporates in mean field the role of a medium that leads to the MSW effect. The 
neutrinos will interact via the weak interactions with electrons produced via evaporation 
of microscopic black holes. Any rigorous discussion of such a process would involve a full 
theory of QG which is not available currently. In the next section some semi-classical ar- 
guments from black hole physics are summarised which motivate this possibility. Of course 
for such a medium it is also necessary to incorporate fluctuations and this will be investi- 
gated at length in the next section through the introduction of a A2 with a specific double 
commutator structure. 

As in the scalar case only weak fluctuations h^^'^ around the flat metric r/^'^ are considered 
and as for that case we will consider the form of g^^ in ()2.4|) . The lagrangian C f for a Dirac 
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particle of mass ruf (in standard notation) is (see, for example, 28|) 



(3.1) 



where h = h^^rj^y (= a\ — a\ — a\ + a\ + 2(ai + 04)). The total lagrangian will have 
contributions from electron and muon neutrino spinor fields \l/e and in the form of ()3.1|) 
together with a Dirac mass mixing term (proportional to iTLf.^) and a MSW interaction. On 
writing 



X 



(3.2) 



where x ^^id (f) represent Weyl spinors, our total Lagrangian, including the mixing and MSW 
terms, becomes j29| 

£ = (1 + hi) [xlidoXe + xUlidlXe + ^l^^O^e - 0lo-l^<9l0e) 

- 2 - hai)doXe + xK^sl - h2cri)diXe) 

- (0t(6,i + 63ai)9o0e + 0l(&3l + &2ai)9i0e) + {e ^ /i} (3.3) 

- (1 + \h){me^,{xl<P^Ji + ^J^Xe + xj^e + ^IXm) + '^0l0e) 

- (1 + hi)m^{xl<Pe + 0lXe) + {e /i} 

Here V is the coupling which represents an MSW effect and is proportional to the density of 
the microscopic black hole density. Moreover, for convenience, we have made the definitions 



hi = a\ + 2ai — 



hi = an 



2a4 



(3.4) 



63 = aitts + as - 0204 - a2- 



29| but now in the presence of a stochastic 



We follow the basic procedure presented in 
gravitational background. In the absence of V the mixing matrix U has the same form as 
in the last section with 

tan (2^) = ^^^ (3.5) 

and so 

^2 




COS 9 sin 9 
— sin 9 cos 9 




(3.6) 
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and 




cos 6 sin 6 
— sin 6 cos 




(3.7) 



This results in 



C = il + ^h){x\iido + iaiidi)xi + X2i^do + aiidi)x2 
+ (p\{ido - aiidi)4>i + (pli^do - criidi)4>2 

- miixlh + - m2{xl<p2 + 02X2) 

- V{cose(pl + sm9(f)l){cose(pi + sm9(l)2)) (3.8) 

- ^(Xl(^ll - &3Cri)5oXl + Xlih^ - &30-l)<9oX2 

+ xlihl - h2a{)dixi + xk^al - h2a{)dix'2) 

- ^(0l(&il + hai)dQ<Pi + 05(6il + hai)dQ(t)2 

+ 0l(63l - &2ai)5l0l) + 4(^3l - &2ai)ai02). 

Owing to translation invariance for the MSW medium in mean field V is constant and we 
make an expansion of the fields in terms of helicity eigenstates 

0, = J]e^^-^{(P:(A;,t) + iV^(A;,t))a(fc) + (P^(A;,t) + iV^(fc,t))/3(fc)} 

k 

X, = Y,(^'''''{{Qiik,t) + Mi{k,t))a{k) + {Q^p{k,t) + M;{k,t))Pik)} (3.9) 

k 

where the motion is in the x-direction, P^, (with fi = a,(3) are positive frequency and 
iV^, are negative frequency field components. The properties of the helicity eigenstates 
can be summarised by the relations 

aik/3{k) = -k(3{k) ail3{k) = -(3{k) (3.10) 
aika{k) = ka{k) =^ aia{k) = a{k). 

On substituting the expansions (|3.9|) into the equations of motion ()B2|) and taking the 
projection of the equations of motion onto positive frequency and negative helicity states 
we obtain 

(1 + i/i) {{ido - k - Vcos^ 6) P^{k,t) - niiQl{k,t) - V cos 9 sin 9P^{k,t)) 
ih - h) Pl{Kt) + ^(63 - h2)Pl{Kt) = 

12 



(1 + ^h) (iQlik, t) + kQlik, t) - miPl{k, t) 

-^ih + h)Q'p{k, t) + ^(63 + b2)Q'^{k, t) = 
(1 + ^h) {{ido -k- V sin'^ 9) P^{k,t) - m2Ql{k,t) - V cos 9 sin eP^{k,t)) 

i{b^-bs)Pi{k,t) + ^ 

(1 + h) {{tdo + k)Ql{k,t) - m2Pl{k,t)) 



~ih - h)Pl{Kt) + -(63 - h2)Pj{Kt) = (3.11) 



2 

+ h)Ql{k, t) + ^(63 + h)Ql{K t) = 

We seek solutions with time dependence e^'^*. This leads to an eigenvalue equation for E (cf 
Appendix B for details). As with the scalar case, to find the flavour oscillation probabihty 
it is necessary to compute (e*^'^^"'^^^*). Gaussian integration gives 

/ 2 u-B-i-u 

^4^g-a.B.a+«.a ^ (3.12) 
Vd^ 

where, in our case, 

u = (i^i!^i-^^i + i2Wcos2^,z^^:^i-^^i + zl^icos2^ (3.13) 
2k 2k 

_.K^, - iVtcos29, iK^t) (3.14) 
2k 2k ' ^ ' 

and the components of the symmetric matrix B are 



= l_,,f K-^i) _4yfeeos2g 

B,2 = ^21= it (^^^^8^-^ cos 2^), 

B,, = ^31 = it (^^^^^^^ + ^cos2^ 

= = {^^^^j^ + Vcos29^ , 

= - + ^f^^^i^ + ycos2^y (3.15) 
(72 2 \ k J 

S23 = ^32 = - T^cos2^ - -^—^ 1 , 



-B24 — -B42 — 



2 V 2k 

it{m\ — m\) 
8^ ' 



13 



= —-^tVcos29, 
0-3 2 

it ( m? — 

S34 = S43 = — - + Vco^2Q 



4:k 



-B44 = — • 

(74 



These expressions have been obtained in the physically relevant limit k'^ ^ m\,m\ and 
|T| ^ 1 where T = i. ■ On using these relations and substituting into eqn. ()3.12|) we 



find 



|/g«{'^l-'^2)i^ 



_ 1 (^-iait(^ '"'i;'"'' +V COS 2e) +1^ ,os 2^) -if^V cos 29 

(m?-mi^)^ ,„ , 2V cos 2e{m?-mT) „ „ „ 9 
1 2fc^ (9<Tl+<72+<73+<74)+ P a^(12<7i+2<72-2<73)h2 

xe V / (3.16) 



where 



-2, 



0-2 



i (m2 + T(l + cos2^)(m2-m2) + T2(m2-m2)sin2 2^) 
I (m^ + T(l - cos2^)(m2 - ml) - T^{m\ - ml) sin^ 26) 



(3.17) 



There is again a suppression of the oscillations which is gaussian with time and also the 
oscillation period is modified in an interesting way which depends both on the square of the 
mass differences, the mean density of microscopic black holes and the effects of back-reaction 
on the gravitational metric. 

Although not done explicitly here, the analysis of the effect of stochastic quantum fluctu- 
ations of the background space-time for the case of Majorana fermions leads to qualitatively 
similar results. 



IV. SPACE-TIME FOAM MODELLED AFTER THE MSW EFFECT 
A. MSW-like effects of stochastic space-time foam medium 

In the suggestion that the observed mass differences between neutrinos are generated 
by a sort of stochastic space-time foam has been proposed. If microscopic charged virtual 
black/white hole pairs were created out of the vacuum then information loss would be 
induced and the subsequent Hawking radiation would produce a medium with stochastically 
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fluctuating electric charges. This radiation would have a preponderence of electron/positron 
pairs (e e) (over other charged particles (muons, etc) from kinematics) and the 'evaporating' 
white hole could then absorb, say, the positrons. According to the Standard Model of particle 
physics, the resultant electric current fluctuations would interact more strongly with z/g 
rather than z/^, and lead to flavour oscillations, and hence, effective mass differences, for the 




neutrinos. This parallels the celebrated MSW effect [22, |23[ for neutrinos in ordinary media. 

From semi-classical calculations there is a significant difference between neutral and 
charged black holes. As neutral black holes evaporate they become less massive and there is 
an increase in the rate of evaporation. Consequently they have a short lifetime. The force 
on a neutrino z/ due to the emitted electron-positron pair is [2^ J2a Gav^a where rio- is the 
particle density of species a in the medium and 

G^^ = % - S^^) {5,,^ - 5,^) (1 + 4sin2 Ow)] + O (4.1) 
V2 \f^w / 

and my/ is the mass of the charged weak boson and 9w is the weak angle. If = n-e then the 

force on a z/g would vanish to O ( ) . Similar subdominant terms are produced for other 

flavours of neutrinos and so neutral black holes would have an equivalent interaction with all 

flavours of neutrinos. On the other hand charged (Reissner- Nordstrom) black holes of charge 

and mass 971 emit electron-positron pairs for 971 > 15 but as 9Jt ^ 13, the extremal black 

hole limit, the surface gravity k — > and evaporation ceases ( see e.g. |3l| and references 

therein) . 

The limiting behaviour of near extremal charged black holes can be made more precise 



from field theoretic studies of black holes 



by actually bounding the number A^^^^ of 



massless (scalar) particles (or pairs of particles/antiparticles) created in a state represented 
by a wavepacket centered around an energy ujq: 



2 c{ujo)\t icoo)\^ 
(2n7r) 



Here c{ujo) is a positive function, > is an arbitrary but large power, £, m are orbital 
angular momentum quantum numbers (arising from spherical harmonics in the wavefunction 
of the packet), and 2mi, n being a positive integer, is a special representation of the retarded 
time in Kruskal coordinates The wavepacket has a spread e in frequencies around ujq^ 
and in fact it is the use of such wavepackets that allows for a consistent calculation of the 
particle creation in the extremal black-hole case. From the expression ()4.2|) . we observe that 
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since 2mr represents time, the rate of particle creation would drop to zero faster than any 
(positive) power of time at late times. The limit of extremality is obtained by means of 
certain analyticity properties of the particle creation number j3l|. In the expression (|4.2p 
t(co'o) denotes the transmission amplitude describing the fraction of the wave that enters the 
collapsing body, whose collapse produced the extreme black hole in [s^. 

In the case of space-time foam, we have currently no way of understanding the sponta- 
neous formation of such black holes from the QG vacuum, and hence in our case, it is an 
assumption that the above results can be extrapolated to this case. In such a situation, then, 
t(c<Jo) would be a family of parameters describing the space-time foam medium. From the 
smooth connection of non-extremal black holes to the extremal ones, encountered in string 
theory js^, we can also conclude that near extremal black holes would be characterised by 
relatively small particle creation rate, as compared with their neutral counterparts. Hence 
black holes which are close to being extremal have long lifetimes. Furthermore when a 
charged black and white hole pair is produced, the absorption of the positron by the white 
hole leaves electrons to preferentially interact with the electron neutrinos. Hence the flavour- 
favouring medium is characterised by charged black/white hole configurations. This flavour 
bias of the foam medium, which could then be viewed as the "quantum-gravitational ana- 
logue" of the MSW effect in ordinary media. In this sense, the QG medium would be 
responsible for generating effective neutrino mass differences [l7|. Since the charged-black 
holes lead to a stochastically fluctuating medium, we shall consider the formalism for the 
MSW effect in stochastically fluctuating media j2J], where the density of electrons replaces 
the density of charged black hole/anti-black hole pairs. It should be stressed, however, that 
we have no way of rigorously checking the required extrapolation to microscopic black holes, 
with the present understanding of QG. However, we shall argue later in this paper, one can 
already place stringent bounds on the portion of the neutrino mass differences that may be 
due to QG foam, as a result of current neutrino data. 
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B. Two Generations of Neutrinos 

Following the MSW formalism, it was proposed in that the stochastically fluctuating 
media caused by the space-time foam can give a mass square difference of the form: 



where k is the neutrino momentum scale and (n^i^ir)) is the average number of virtual 
particles emitted from the foam. These flavour violating effects would contribute to the 
decoherence through quantum fluctuations of the foam-medium density by means of in- 
duced non-Hamiltonian terms in the density matrix time evolution. In this paper we model 
this foam/neutrino interaction by analogy to the MSW interaction Hamiltonian and follow 
corresponding procedures to calculate the relevant transition probabilities. Moreover, QG 
induced Gaussian fluctuations of energy and oscillations lengths may be distinguished from 
the corresponding ones due to the conventional uncertainties by their energy dependence: 
conventional effects decrease with increasing (neutrino) energy, whilst QG effects have ex- 
actly the opposite effects, increasing with energy. 

In keeping with our analysis of the effects of Ai, and for simplicity, we restrict ourselves 
to the case of two generations of neutrinos which suffices for a demonstration of the generic 
properties of decoherence. We take the effective Hamiltonian to be of the form 



where Hj is a 2x2 matrix whose entries depend on the interaction of the foam and neutrinos 
and H is the free Hamiltonian. For the purposes of this paper we take this matrix to be 
diagonal in flavour space. Although we leave the entries as general constants, a^- , we expect 
them to be of the form cx G'jv?7,^j^(r); so we write Hj as 



where the foam medium is assumed to be described by Gaussian random variables [16|. We 
take the average number_of foam particles, {n^i^it)) = riQ (a constant), and {nlf^{t)nlf^{t')) ~ 
^fnl^it-t'). Following |2J] we can deduce the modifled time evolution of the density matrix 

as 



(4.3) 




(4.4) 



d_ 

dt 



(P) 



(4.5) 
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where (...) represents the average over the random variables of the foam. The double com- 
mutator is the CPT violating term since although it is CP symmetric it induces time- 
irreversibility. It is also important to note that A2 here is of the Markovian-Liouville- 
Lindblad form for a self-adjoint operator. This is as an appropriate form for decoherence for 
environments about which we have little a priori knowledge. In the CPT violating term we 
can require the density fluctuation parameter to be different for the anti-particle sector from 
that for the particle sector, i.e. Cl ^ Q, while keeping = '^o the same in both sectors. 

Physically this means that neutrinos and antineutrinos with the same momenta, and hence 
interacting with the same amount of foam particles on average, will evolve differently; this 
is a result of CPT violation. 

We expand the Hamiltonian and the density operator in terms of the Pauli spin matrices 

(with Y = I2 the 2x2 identity matrix) as follows 

Heff = i2(^, + noh'f-^, P = i2p^Y- (4-6) 

(where H^ff = H + uqHi). We find that 

m| + m| mj-ml 
K = ^^^MO + ^^^.3 (4.7) 



and 



iT-oh'f, = t, ^/^o + (^i'e - (^u^) sin 26* + {a^^ - a^^) cos 26* ^,,3. (4.8) 



The master equation in (j4.5p simplifies to 

3 

Pi = J2^ijPj- (4-9) 

for Z = 1, . . . , 3 (see Appendix C for further details). The pure state representing i/g is given 
by 

(p)(-) = + sin {26) I + cos (2^) | (4.10) 
and the corresponding state for is 

{Pf^^ = il2 - sin (2^^) I - cos (2^^) |. (4.11) 

If (p) (0) = (p)^'^"^ then the probability P^^^y^ {t) of the transition z/g is given by 

P..^.At)=Tr{{p){t){p)^''^^). (4.12) 
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In order to study decoherence we will calculate the eigenvectors and corresponding 
eigenvalues Aj of C to leading order in In terms of auxiliary variables U and W where 



U = (a^, - a^J cos {26) + 



mr — mo 



2k 



(4.13) 



and 



W = (a^, - a^J sin (2^) 



it is straightforward to show that 



^1) 



_u_ _.VWTW 
'W' vv ' 



(4.14) 



(4.15) 



and 



Ai ~ -fi^ (Wcos(2^) -Wsin(2^))% 
^ n2 

A2 ~ -zVW2 + W2 - — (W^ + + (W COS (20) + W sin (20))^) 



A3 ~ i VW^TW - — ^ + (W COS (20) + W sin (20))^) . 
In fl4.9p the vector "p* (0) can be decomposed as 



(4.16) 



(0) = biT'^) + bsT'^) + ba"? 



(3) 



(4.17) 



with 



and 



Hence 



and so 



bi 



cos (29) +UW sin (29) 

cos(20) -ZYWsin (20) 
2 (W2 + w2) ' 



1 



p (t) = - (bie^i*-?'^) + baT'^) + bsTt^) + I2) 



P. 



it) 



1 - sin (20) jbiel'^e^i* + ba (ef ^e^^* + eS'^e^^*) } 
-cos (20) {bie^'^e^i* + bs (cf e^^* + ef e^^*)} 



(4.18) 



(4.19) 



(4.20) 
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On writing A = Oj,^ — and 5k 



it) 



^2]^^ 1 Pfe^f^, (t) readily simplifies to give 
Ti (t) + T2 (t) 



2 (A2 + 52 + 24 A cos (20)) 



(4.21) 



where 



ri(t) = (A + cos (20) 6kf (l-e- 



sin^(2e)5lt 



(4.22) 



and 



r2(t) 

= 6lsm^ (20) < 



cos 



X exp 



^2 



24 A cos (20)t 



(4.23) 



-f (2 (A + cos (20))' + 61 sin' (20)) t 



Since we are concerned with relativistic neutrinos, we have t = x (in natural units) and 
we can use this to put our expression in terms of the oscillation length, L. The exponent in 
the damping factor in ()4.21|) has a generic form 

exponent oc i7'/ (0) L 



with / (0) = {A + Sk cos (20)) +2^1 sin (20) or — — j^-^. Hence the damping is directly pro- 
portional to the stochastic fluctuations in the medium. The limit 6^ characterises the 
situation where the dominant contribution to neutrino mass differences is due to space-time 
foam (Q|. The damping exponent should then be independent of the mixing angle for consis- 
tency. Indeed we find the purely gravitational MSW to give exponentg^avitationaiMsw ^ A^L 
which is independent of 0. However this stochastic gravitational MSW effect, although capa- 
ble of inducing neutrino mass differences, gives an oscillation probability which is suppressed 
by factors proportional to S^. Hence the bulk of the oscillation is due to conventional flavour 
physics. 



C. Comparison with decoherence from conventional sources 

In experiments with neutrino beams there is an uncertainty over the precise energy of 
the beam (and, in some cases, over the oscillation length), which can destroy coherence, as 
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discussed in [25". There are also small effects due to the wavepacket nature of the incoming 
neutrino state. The coherence length associated with the latter is typically much larger 
than L and so a plane- wave approximation is sufficient. Below we first review the situation 
briefly, for the benefit of the inexpert. 

In refs. jssf 3| the following expression for the neutrino transition probability has been 



considered: 



n n / 

P^^p = P^p{L, E) = 5o.p- 4 J] J] ^{U*^,UpaUc.bU;,) sin^ I 

a=l 6=1 ^ 



AE 



a<b 



^=1 1 V / 



a=l 6=1 
a<b 

where L is the neutrino path length, E is the neutrino energy, n is the number of neutrino 
flavours, and Am^^ ( = — ml) and Uaa as before is the mixing matrix. As there are 
uncertainties in the energy and oscillation length, in refs. 
the L/E dependence was taken. This average is defined by 

(P) = / dxP{x)^^e-^^. 
J-oQ crv27r 



26| a gaussian average over 



where x = j^, I = {x) and a = ^/ {{x — (a;))^). Furthermore if L and E are independent then 
/ = {L/E) = {L)/4:{E) (for highly peaked distributions) and one obtains for the averaged 
expression 

n n 

P^0{L,E) = 5^p- 2 ^^5R([/*,C/^,[/,,t/|,)(l-cos(2IAm2,)e-2-'(^™^^)') (4.24) 

a=l 6=1 

a<b 
n n 

- 2 Y.Y.^^^*o.aUpaU^bU^p^,)s\n^ {2l^ml^)e-^''"^''<^^\ a, P = e, ^,t, 

a=l 6=1 
a<b 



It should be noted that / has to do with the sensitivity of the experiment and a the 
damping factor of neutrino oscillation probabilities. A pessimistic (less stringent) and an 
optimistic (more stringent) upper bound for a (obtained from a first order Taylor expansion 
of X around (E) and (L) ) can be given 

• pessimistic: a ~ Ax = A^ < AL||f |^^^^^^^^^^^ + AE||||^^^^^_^^^^) 
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• optimistic: ^<^^{^y+{^y 

For the case of two gene rations, using this procedure, the transition probabihty between 
flavour eigenstates is j25| 



{P..^^J = -sin^2e ^1 - e--^---^^ cos J (4-25) 

Owing to the averaging over Gaussian fluctuations, 14.251 shares one characteristic with 
the back reaction effects of Ai (discussed earlier) viz. the dependence of the decoher- 
ing decay and is dissimilar to the L dependence of the space-time foam (as modelled by 
the gravitational MSW effect). This clearly, in principle, is a way of distinguishing the 
MSW type effect. Although typically experimental data make allowances for systematics, 
it is interesting to consider whether for a given L the magnitude of the decoherence effect 
may be assigned to conventional sources. When one compares the damping factors of the 
conventional averaging and our MSW effect we get 



2a^(Am?2)^ 
which we can express as 



(A + 6k cos (2^))^ + ^61 sin^ (2^) 



(4.26) 



(A + 6k cos (2^))^ + ^61 sin' {29) = ^-^^Lr' (4.27) 
z oh/ 



where i" = ^ + ^ for the pessimistic case or r = y i^)"^ + optimistic 
case. For decoherence due to standard matter effects with L ~ 12000-ft"m, r ~ C>{^), 
E ~ 0(l)GeV , Amj^ ~ O{l0-^)eV' and ~ 1.5 x lO^^GeV"^ one obtains -fatmjake{= 



^-^^Lr') < 10-24 GeV. 



8^2 

It is worth pointing out here that such a small order of magnitude is of a similar order to 
that found ni^quantum gravity decoherence suppressed by a single power of Planck mass 0, 
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34| 



InQ the cases for the decoherence damping factor being of the form 7 = 70 (fj^)"^; 
with 7o as a constant, has been analyzed for the n = 0, —1, 2 cases (a more pessimistic view 
is presented in with 7 = for which there is no experimental sensitivity at least 

in the foreseeable future). An effect of a similarly miniscule order appears to characterise 
also cosmological decoherence, i.e. the decoherence due to theffuture) horizon in de Sitter 



space, in the case of a Universe with a cosmological constant j7l. l36l| ) . 
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In order to investigate experimental signals of quantum gravitational decoherence it will 
be necessary to distinguish genuine quantum gravity effects from the above "fake" ordinary- 
matter effects through the dependence of the respective transition probabilities on the energy 
and oscillation length. Indeed, it is expected, at least intuitively, that the "fuzzyness" of 
space-time caused by quantum-gravity-induced stochastic fluctuations of the metric tensor, 
would lead to effects that are enhanced by the energy of the probe, i.e. the higher the energy 
the greater the back-reaction on the surrounding space-time fluid. Such an expectation is 
confirmed in detailed microscopic models of the so-called D-particle foam jsl. Then, in such 
cases we may write in a generic way 

AL AE E Y , , 

for some positive integer a > 1, and some coefficient, (3. For this case we would have 
r (3 [j^^ then from the gaussian average we would have 



n^(A + 4cos(2^)r + -5^in^(2^)~^^^/3^^^J L (4.29) 

For the specific model of D-particle foam of ref. jj] a = +1, and Mqg ~ Ms/qs with Mg 
the string scale and gs < ^ the (weak) string coupling. 

Since for the oscillation length L, ~ from ()4.29|) and the above analysis, it 

becomes clear that genuine quantum gravity effects in some models are characterized by 
damping factors which are proportional to E"^"", a > 1, and thus are enhanced by the energy 
of the probe, leading to significantly more damped oscillations for high energy probes as 
compared to the low-energy ones. This is to be contrasted with the conventional effects, due 
to the passage of neutrinos through matter, which are diminished with the energy Q]. 
Although in the presence of A2, as shown in |[32], the CPT operator cannot be defined, the 

/ A pCPT \ (decoh) 



CPT violating difference between neutrino and antineutrino sectors |2 



V ^l3a J 



p(decoh) 

^^aecoh) — 1 vanishes unless the decoherence coefficients between particles and antiparticles 
are distinct, considered in Here the superscript decoh denotes the decohering 

piece of the relevant probability. In the case of different decoherence coefficients between 
particle and antiparticle sectors, the QG induced difference AP^^^ would either increase or 
decrease with energy, at least as fast as a Gaussian, depending on the relative magnitudes 
of the decoherence parameters in the neutrino and anti-neutrino sectors. In contrast the 
conventional matter induced CPT difference saturates with increasing E. In this way, at 
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least in principle, the two effects can be disentangled. It must be noted, though that, as 
seen from ()4.29|) the proportionality coefficient /^^(Amfg)^ accompanying {E/MqgY'^{L/ EY 
in the decoherence exponents is very small (for natural values of /3, although in principle 
this is another phenomenological parameter to be constrained by data). Hence, for this 
particular model of QG decoherence, appreciable effects might only be expected in situations 
involving very high energy cosmological neutrinos. In view of this, the analysis of high-energy 



neutrinos performed in 



2l|, which was based only on conventional Lindblad decoherence. 



needs to be repeated in order to incorporate the above effects. 



V. CONCLUSIONS AND OUTLOOK: PRELIMINARY DATA COMPARISON 

It is hoped that decoherence due to quantum gravity can be confirmed or ruled out by 
physical observation. We will make a few remarks concerning possible conclusions from 
data from reactors and the atmosphere. Different approaches have been used in examining 
transitions of atmospheric neutrinos. As mentioned above, more pessimistic expressions for 
damping factors such as 7 = ^EvJp have been presented [35| . However, more optimistic val- 
ues can be obtained. In [3J] a phenomenological analysis is done for the case of atmospheric 
neutrino transitions (z/^ ^ Vr). They obtain upper bounds to the decoherence parameters 
and find that the Super-Kamiokande data can be a be a good probe into quantum gravity 
induced decoherence. They discuss three possible energy dependencies of the decoherence 
parameter, in particular 7 = 70 {E/GeV)^ with n = —1,0,2, with 7n a constant, and the 
subsequent constraints. The controversial data obtained by LSND 38[, if confirmed by fu- 
ture experiment (for instance MiniBOONE), could provide important data which may lead 
to evidence of space-time foam interacting with antineutrinos. 

We would now like to mention briefly some preliminary attempts to constrain the models 
presented here by means of currently available neutrino data. In a recent work, 0| we 
have presented a fit of a three-generation (completely positive) Lindblad decoherence 



model for neutrinos with mixing to all the available data, including the 



;SND result in the 



antineutrino sector. In contrast to the manifestly CPT- violating fit of ^|, which attenapted 
to explain the LSND result from the point of view of CPT-violating decoherence, in it 
was assumed that the decoherence coefficients were the same in the particle and antiparticle 
sectors. The best fit that was obtained showed that only some of the oscillation terms in 
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the three generation probabihty formula had non-trivial damping factors; moreover over an 
oscillation length the exponent of such non-trivial damping, T) ■ L, satisfied 



^ 1.3-10-2 

V = , (5.1) 

in units of 1 /km with L = t the oscillation length. 

In the light of ()5.1|) it is possible to analyse 3] the two types of theoretical models of 
space-time foam discussed in sections III and IV of the present paper. The conclusion is 
that models incorporating stochastically fluctuating MSW-like QG media as in 1)4.211) cannot 
provide the full explanation for the flt. Indeed if the decoherent result of the flt ()5.H) was 
exclusively due to such a model, then the pertinent decoherent coefficient T> in the damping 
exponent, for, say, the KamLand experiment with an L ~ 180 Km, would be \V\ = Q'^A'^ ~ 
2.84-10-21 GeV (note that the mixing angle part does not affect the order of the exponent). 
Smaller values are found for longer L, appropriate to atmospheric neutrino experiments. In 
this context the L independence ofV-L, as required by ()5.1)1 . may be interpreted as follows: 
fl4.21|) suggests that we write A = where ^ <^ 1 parametrises the contributions of the 

foam to the induced neutrino mass differences. Hence, the damping exponent becomes in 
this case E^^VL^ {Ani^Y ■ L / E"^ . Thus, for oscillation lengths L (since ~ Aw? jE) one is left 
with the following estimate for the dimensionless quantity ^"^ Aiv?D? / E ~ 1.3 - 10"^. This 
implies that the quantity V? is proportional to the probe energy E. Since back reaction 
effects, which affect the stochastic fluctuations V? , are expected to increase with probe 
energy E, this is not an unreasonable result in principle. However, due to the smallness of 
the quantity Arn? jE^ for energies of the order of a GeV, Aw? ~ 10"^ eV^ and ^ <^ 1), 
we can conclude that ^ ^ in this case, would be unrealistically large for a quantum-gravity 
effect in the model. We remark at this point that, in such a model, we can in principle 
bound independently the f2 and A parameters by also examining the period of oscillation. 
However in this example, Aaey^ ^ A12 and so the modiflcation in the period is too small to 
be detected. 

The second model ()3.16|) of stochastic space-time can also be confronted with the data. 
In this case ()5.1)) would imply for the pertinent damping exponent 

[w? — mT?, , 2V^cos26'(m? — m^) , A , 

^ ' ' ■(9(7i + a2 + (T3 + a4) + ^(12ai + 2^2 - 2^3) 



1.3-10-2. (5.2) 
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Ignoring, for simplicity, subleading MSW effects from V, and considering oscillation lengths 
t = L ^^-^^^i^ , we observe that the experimental fit (jS.ll) . may be interpreted, in this case, 
as bounding the stochastic fluctuations of the metric (|2.4p to 9cri + (T2 + (13 + (14 ~ 1.3. ■ 10^^. 
Again, this is too large to be a quantum gravity effect, which means that in this model the 

contributions to the damping, (j3.16p . due to stochastic fluctuations of the space-time 
metric cannot be the sole explanation of the fit of js^. 

The analysis of 3^ also demonstrated that, at least as far as the order of magnitude of 
the effect in fj5.1|) is concerned, a reasonable explanation is provided by Gaussian-type energy 
fluctuations, due to standard physics effects, leading to decoherence-like damping of oscilla- 
tion probabilities of the form (j4.25j) . The order of magnitude of these fluctuations, consistent 
with the independence of the damping exponent on the oscillation length L (irrespective of 
the power of L), is 

—^-1.6-10"^ (5.3) 
E 

if one assumes that this is the principal reason for the result of the fit. 

However, not even this can be the end of the story, given that the result fl5.1|) applies 
only to some but not all of the oscillation terms; this would not be the case expected for 
standard physics uncertainties (j4.25|) . The fact that the best fit model includes terms which 
are not suppressed at all calls for a more radical explanation, and so the issue is still wide 
open. It is interesting, however, that the current neutrino data can already impose stringent 
constraints on quantum gravity models, and exclude some of them from being the exclusive 
source of decoherence, as we have discussed above. 

We reiterate that, within the classes of stochastic models discussed, one can safely con- 
clude space-time foam can be at most responsible only for a small part of the observed neu- 
trino mass difference, and certainly the foam-induced decoherence cannot be the primary 
reason for the result of the best fit ()5.H) . obtained from a global analysis of the currently 
available neutrino data. Of course, it is not possible to exclude other classes of theoretical 
models of quantum gravity, which could escape these constraints. At present, however, we 
are not aware of any such theory. 

In the near future we plan to make a more complete and systematic comparison of our 
new formulae, especially those derived in sections II and III, with all experimental data 
available and perhaps arrive at new constraints. 
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APPENDIX A: SCALAR PARTICLE AVERAGES 

For integration over metric fluctuations we shall use the formula 



(Here the a's are assumed to be in the range (—00, 00) and the form of B must be such that 
'convergence' of the integral is assured.) 




B = 



H - it(f(mi) - f (777,2)) 



For simplicity we define 



T = f(777i) - f(7772) 



and 




c 



777 



2(A;2 + m2)-3/2_^2^A;2 + m2)-3/2 



So we can write 



^11 




ml + 2k'^ ml + 2k'^ 



•^23 — 




and the remaining J^ij — 0. 
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Putting this information together we find 



( ^ 



iht 








2d ' 
-ikiht 
2d 







-ikHt 
2d 
J_ 








^ - hk'^ct / 

(74 Z / 



Ui = 
Ui — 

i.e. u — 



-itb 

d 

ith f -1,0,0,-- J 



where 



So we obtain 



detB 



1 



16(71(72(73(740?^ 



^PxP2 



Pi 



P2 = 



+ 2idLa2kH - 2iba2dH + h'^kH'^a2(Tz 
Ad^ + 2idha2{k^ - d)t + h^k'^a2(T^t^ 

- 2id^{k^5a4 + 2bai)t + bk^aia^ibk^ - 2d^c)t^ 

1 



detH 



det H 



1/2 



16rf 



(71(72(73(74 
1/2 



det By \P1P2; {P1P2Y''' 



Now 



_ 2aibt{k^a4{bk^ - d^d)t - 2id^) 

~ AdP + baik^a^ibk^ - 2id?)t'^ - 2i{k^aTi + 2ba^)d? 
""2 = 0, f 3 = 

-2{bait + 2i)bda4kH 



Ad^ - 2zd'^{k^ca4 + 26cTi)t + bk^aia^f^ibk^ - 2cd^) 



exp('ii • v) — exp ( — 
X2 
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where 



APPENDIX B: DIRAC PARTICLE AVERAGES 



The equations of motion which follow from (j3.9p are 

(1 + -h){idQ(j)i — {iaidi(j)i + rriiXi)) — V cos 6 [cos 6 (pi + sin 6*02) 
+ 63^i)9o0i + (63I + h2(Ji)di<Pi) = 

(1 + ^h){idoXi + iodiXi - rriicpi) 

- 630-i)<9oXi + - &20-i)<9iXi = 

(1 + ^/i)(i9o02 - ^cri<9i02 - "^2X2 - ^ sin 6^ (cos 6^01 + sin6'02)) 

-^((62! + &3fXi)9o02 + (&3I + &2C^l)5l02) = 
(1 + ^/l)(i9oX2 + iOidxX2 - "^202) 

--((62I - &3f^l)5oX2 + - &20-l)5lX2) = 

On using ()3.9p in (|B2|) we have 

Q\{k,E) 

\Q}ik,E)J 
where M is a 4 x 4 matrix with components 



M 



M12 

Mi3 = - 

Mi4 = 



E(l + i/^-i(6i-63; 



1 



1 



(1 + -h)k - (1 + -h)Vcos'{e) 
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-(1 + -h)mi 
-1/(1 + i/i) sin(^) cos(^) 



.{l + -h)rm 
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M23 = M24 = 

M31 = -(1 + ^/i)V^cos(^)sin(^) 

M32 = 

M33 = E(^l + h- ^(61 - 63)) - A; (^1 + l/i + 1(62 - 63)) - (1 + h)Vsm\e) 

M34 = -(1 + l/i)m2 

M41 = M42 = 

M43 = -m2(l + l/i) 

M44 = E(^l + h- l(6i + 63)^ +k(^l + h + 1(62 + 63; 

~ 1 2 

Using these equations one can eliminate by substitution to obtain 

AT I I = (B3) 



where 



Ml = Mu + ^rm{l + h) 
IVI 22 ^ 

A/12 = —V sin 6* COS 

N21 = M31 (B4) 

U22 = M33 ' ^ 

M44 

2 

We take the momentum k to be very large, and so we write ~ /c + We make the 
substitution 

= zo + ^ Zitti + ^ Zijaittj (B5) 

i ij 

and expand the components of N in terms of the stochastic parameters Oj. This allows us 
to use the condition det = to find the Zi terms. There are two solutions of labelled 
by z^ and zf. 

We use (jX)) to evaluate 

/det 
rf^aexp(-a-S-a)e'("i-"^)*-^ (B6) 
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with 



and 



it 



u 



2k ' 



^4) 



B 



it ( 



'11 



12) 



' 2k (>^13 ^ 



13J 



^12 J 



CT2 ^(^22 ^22) A: 



2k 

i|(^23 ~ ^23) 



'2^1^14 



'147 



2fcV^24 ^' 



24 J 



■2fcV^13 %3J 



■2fcV^23 ^23 J 



^33)1 2fc(^34 ^34) 



,(B7) 



\ ~i|(^14~^14) ~i|(^24 "~ ^24) ~i|(^34 " ^24) ~; ~ "^(^44 ~ ^44)1 / 

On substituting the detailed expressions for Zq and zf it is straightforward to obtain the 
forms in (Hm) and (ITT!)!) . 



APPENDIX C: LINDBLAD DECOHERENCE 



A useful generic form of the Lindblad master equation for a x iV density matrix p is 



where ^ ^ 



kl=l 



The complex N x N matrices Fi ^= F/^ , / = 1, . . . , A^ — 1, together with the identity 
matrix Iat (= Fq) form a basis for a space of complex A x A matrices and so any operator 
D can be written as D = X]^=o ^ i^fiF^^. If {cfcz} is a non-negative matrix, Tr (Fi) = 0, and 
Tr (FiFj) = ^Sij, then the density matrix p evolves in the space of physical density matrices 
[1^ and so probabilities are non-negative. On writing H = ^^=0 ^m-^m have 



2 

j,fc=i fc,/=i 



where 



1 / [F,, [p, FiW + {F,, [p, F,]} + [[F,, p] , F,] 

2 \ +{[Ffc,p],F,} + 2{p,[Ffc,F;]} 
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For N = 2, Fj = Y (where Sj are the Pauh matrices) Oq = ^Tr (O) and Oj = Tr (Osj). 
The master equation of ( 14. 5|) becomes 
d 



dt 



ip) = -i[H + n,Hi, (p)] + 9.'ni {[Hj {p) , Hj] + [Hj, {p) Hj]) 



on noting that 

[Hj, [Hj, (p)]] = -{[Hj (p) , Hj] + [Hj, (p) H j]) 

The non-zero elements of the associated c matrix for ()C5|) are 

2„- 2- 



(C5) 



(C6) 



cii = 2Vt^ [a^^ - ay^) sin^ 26^, 
Ci3 = C31 = 2f2^ (oj,^ — a,^^)^ sin 2^ cos 2^, 
C33 = 2fi^ (ttj,^ - a,,^)^cos^26'. 



(C7) 



On using (j4.6p 



Also 



1 ^ s 



(a 



(C9) 



j,r=l 



Po is independent of time from the structure of ()C5|) whereas Pq (g = 1, 2, 3) satisfies 

d ^ 

~ 2 



^ Cpi {2Sjg5pi - SjpSgi - pj. 



(CIO) 



Using this it is straightforward to show that the C corresponding to ()4.9j) is 
^ A2 cos2 (2^) -U gjj^ (20) COS (2^) ^ 

A2 sin (20) cos (20) W -Q"^ A"^ sin"^ {26) J 
where U and W are defined in ()4.13|1 and ()4.14|1 . 



v 
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